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conduction  by  electron  tunneling  through- thin-film  aluminum  oxide  junctions.  Details  are  presented 
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GENERATION  OF  INTERMODULATION  BY  ELECTRON 
TUNNELING  THROUGH  ALUMINUM  OXIDE  FILMS 


INTRODUCTION 

The  Navy’s  Fleet  Satellite  Communication  System  (FLTSATCOM)  consists  basically 
of  three  synchronousorbit  satellites  with  UHF  links  to  various  earth-surface  terminals. 
Each  satellite  and  surface  terminal  simultaneously  operates  both  transmitter  and  receiver 
units  in  close  physical  proximity.  The  Navy’s  transmit  and  receive  bands  are  approxi- 
mately 240  to  270  MHz  and  290  to  320  MHz  respectively.  These  two  carrier  bands  are 
separated  by  only  20  MHz.  For  such  a system  to  operate  successfully,  any  undesirable 
frequency  mixing  caused  by  nonlinear  components  (intermodulation  (IM)  generation) 
must  be  minimized. 

The  severe  magnitude  of  this  problem  can  best  be  illustrated  by  calculating  the 
acceptable  IM  power  level  at  a receiver  for  typical  operating  conditions;  for  example,  for 
a satellite  transmitter  power  of  ^ 100  W (+50  dBm)  the  space  attenuation  alone  is  ^ 175 
dB  (-  125  dBm,  or  down  to  3 X 10~^^  W)  for  either  the  uplink  or  downlink  signals.  At 
the  satellite  receiver  the  uplink  signals  should  normally  be  ^ 15  dB  above  any  interference 
or  noise  level.  This  requires  the  interference  level  to  be  below  -140  dBm,  or  lO'^"^  W. 
Hius  the  acceptable  IM  level  at  the  satellite  receiver  is  required  to  be  190  dB  below  the 
local  transmitted  power  level.  This  is  a ratio  of  IM  power  to  signal  power  of  10'^^.  The 
same  requirement  must  also  be  met  at  the  earth-based  terminals. 

This  problem  is  diagramed  in  Fig.  1 by  an  illustration  from  Young  [1].  Here  the 
IM  products  are  presented  in  a simplified  two-signal  analysis.  Consider  two  frequencies 
fi  and  f2  in  the  transmit  band  at  voltage  levels  of  Vi  and  ^2*  input  voltage  Yin 
across  some  network  can  be  written 

Vjn  = Vi  cos  + V2  cos  cj2f-  (1) 

For  a linear  network  the  output  current  will  be 

/-GoVin. 

For  a nonlinear  network  the  output  current  can  be  expressed  as  the  power  series 

/ = GoVi„+GiV«  +G2V3  +....  (2) 
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(a)  Basic  satellite  communication  system 


(b)  SimpUfled  two-signal  analysis 

Fig.  1 — Basic  satellite  communication  system  showing  simplined  two-signal  analysis  for  a 
linear  and  a nonlinear  network.  The  intermodulation-frequency  terms  result  from  a power 
series  representation  of  the  nonlinear  signal-voltage  output  and  constitute  interference  in  the 
receive  bands. 
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The  separate  IM  term  in  the  nonlinear  output  can  be  seen  by  substituting  the  input 
voltage  of  Eq.  (1)  into  Eq.  (2)  to  give 


/=  Go(^l  cosojjf  + V2  COSCJ2O  + 
G 


^1*^2  1 1 
— * ~2  2wit  * ~2  ^2  2002^ 

+ V■^^  V2  cos  (cjj  - CO2)  + V2  cos  (cj^  + C02)j 

'1  1 

— Vj  cos  Sco^t  + Vg  cos  3cj2t  + 3Vi  cos  tjjt 

fvl  vl\  3 

* ^ "4*/  '^l  (2^2  “ 

3 3 

+ 7 Vj  V|  cos  (2aj2  + Wj ) + 7 V,  cos  (2cj^  - o;2) 


Q 

+ 7-  Vj  Vg  (2wj  + 10 2) 


(3) 


The  frequencies  ==  250  MHz  and  ^2  = 270  MHz  have  been  chosen  for  most  of  the 
device  testing  at  the  NRL  IM  test  facility. 

At  these  two  typical  satellite  transmitter  frequencies,  most  of  the  IM  terms  shown 
in  Eq.  (3)  include  frequencies  that  lie  outside  the  satellite  receive  band.  However,  two 
of  the  largest  magnitude  IM  terms  of  frequencies  2/^2  - fi  and  8/2  - 2fy  (290  MHz  and 
310  MHz)  do  fall  within  the  local  receive  band.  (The  source  of  the  latter  IM  frequency 
is  the  term  and  is  not  shown  in  Eq.  (3).)  This  example  illustrates  the  problem 
resulting  from  a narrow  separation  of  the  transmit-receive  bands. 

The  NRL  investigation  is  concerned  with  the  generation  of  intermodulation  due  to 
various  nonlinear  mechanisms  in  normally  passive  hardware  such  as  contacts,  flanges, 
connectors,  and  other  current-carrying  me^  structures  attendant  to  the  transmitter- 
receiver  system.  The  investigation  described  in  this  report  is  speciRcally  concerned  with 
the  generation  of  intermodulation  due  to  the  nonlinear  conduction  by  electrons  tunnel- 
ing through  thin  oxide  films,  particulariy  AI2O3  on  aluminum  structures. 

The  detailed  mechanisms  of  surface-to-surface  contacts  is  extremely  involved  and 
complex  and  remains  incompletely  understood.  Ihere  are  many  interacting  variables, 
such  as  the  density  of  microscopic  contact  points,  number  of  contacts,  contact  pressure, 
oxide  growth,  oxide  fracture  and  regrowth,  metal-to-metal  bonding,  diffusion,  and  time 
and  temperature  effects. 

Despite  the  lack  of  a complete  microscopic  theory,  contact  models  such  as  those 
described  in  the  Philco-Fotd  report  [2]  do  aUow  predictions  adequate  for  assessment  and 
control  of  the  intermodulation  generate  in  real  hardware  contacts.  Because  of  the 
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numerous  contact  mechanisms  which  have  been  identified,  an  experimental  study  of 
tunneling  in  real  contacts  is  not  amenable  to  reproducible  results  and  to  correlation  with 
theory  [3],  Consequently  the  intent  of  this  experimental  investigation  is  to  fabricate 
well-characterized  tunneling  junctions  of  AI-AI2O3-AI,  to  directly  measure  the  inter- 
modulation generated,  and  to  attempt  to  correlate  the  intermodulation  with  the  device 
circuit  parameters  and  with  tunneling  theory.  In  addition  some  preliminary  efforts  were 
made  to  improve  the  conduction  characteristics  of  such  junctions  by  ion  implantation. 


TUNNELING  THEORY 

The  theory  of  electron  tunneling  through  an  insulating  layer  dates  back  to  the  1920’s. 
However  controversy  still  continues  as  to  the  correct  form  of  the  tunneling  equations 
for  thin  films  5 nm).  Some  of  the  principal  difficulties  include  the  following: 

• Applicability  of  macroscopic  parameters,  such  as  the  dielectric  constant,  to  a few 
atomic  layers, 

• Complexity  of  integrals  which  can  presently  be  solved  only  by  approximations 
which  affect  the  accuracy  of  the  final  results  to  an  undetermined  degree,  and 

• Accurate  evaluation  of  contaminant  effects  on  the  electron  scattering  surfaces. 

This  last  effect  has  been  developed  into  a valuable  qualitative  method  for  investigating 
the  vibrational  modes  of  the  contaminants,  but  quantitatively  the  effect  is  poorly  under- 
stood. Thus  comparison  of  the  experimental  results  with  the  electron  tunneling  theory 
is  only  approximately  quantitative. 

The  theoiy  of  electron  tunneling  is  based  on  the  fact  that  the  wave  function  of  a 
free  electron  \l/g  extends  only  a few  nanometers  into  an  insulator  at  any  metal-insulator 
interface.  In  large-dimension  insulators  the  wave  function  quickly  decreases  to  a vanish- 
ingly small  value  such  that  the  effect  goes  unnoticed;  but  if  the  insulator  is  only  3 nm 
(a  typical  aluminum  oxide  film),  the  wave  function  remains  finite  throughout  the 
insulator  and  yields  a measurable  probability  for  the  electron  to  pass  from  one  side  of 
the  insulator  to  the  other.  This  effect  is  shown  in  Fig.  2,  where  the  extended  wave  func- 
tion \pg  on  the  left  is  matched  in  magnitude  and  slope  to  the  decaying  exponential  within 
the  insulator,  which,  before  it  vanishes,  similarly  matches  with  the  extended  wave  function 
in  the  metal  on  the  right. 

Since  temperature  has  a small  effect,  one  needs  only  the  zero-temperature  expression 
for  the  tunneling  current  density  J in  one  direction  [4] : 

The  parameters  are  shown  in  Pig.  2,  where  V is  the  voltage  applied  across  the  junc- 
tion and  is  the  Fermi  level  of  the  electron  distribution  in  the  metal  on  the  left.  P{E) 
is  the  transmission  probability  by  tunneling  for  electrons  of  transverse  energy  E, 


I (Efj  - E)P(E)  dE 


4 


NRL  REPORT  8170 


Fig.  2 — Electron  tunneling  through  the  potential  bar- 
rier of  a thin  insulating  film  between  two  conductors, 
llie  conductors  of  the  sandwich  structure  are  main- 
tained at  a potential  difference  of  V such  that  there  is 
a finite  electron  probability  function  everywhere  to 
the  right  of  the  flrst  interface.  The  resulting  current  is 
a nonlinear  function  of  voltage. 


The  usual  method  for  calculating  the  transmission  probability  involves  the  WKB 
approximation,  which  yields  a transmission  probability  given  by 


?(£■)  = exp 


In  this  expression,  ^(ac)  is  ttie  local  barrier  potential  acting  on  the  electron  within  the 
classically  forbidden  region  of  the  aluminum  oxide,  and  the  integral  is  carried  out  over 
the  width  of  the  forbidden  region  S.  For  the  simple  trapezoidal  barrier  shown  in  Fig.  2, 
the  transmission  probability  from  left  to  ri^t  is  given  exactly  by 


P(E)  = exp|-  -I  Y [eV  + A j + <^0  + - (^  + 00 ■ ■ 


This  exponential  expression  is  sufficiently  complex  so  that  the  last  integral  in  the  initial 
equation  for  J cannot  be  carried  out  exactly.  Forlani  and  Minnaja  [5]  have  expanded 
the  argument  of  the  exponential  about  the  average  of  the  Fermi  levels  and  obtained  an 
approximate  analytical  expression  for  the  tunneling  current  density  J.  As  part  of  this 
investigation  a Taylor  expansion  of  the  Forlani-Minnaja  equation  to  the  fourth  order  in 
voltage  has  been  published  [6] : 


. me  -M 
e 

2x2  h® 


480 


(4) 
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where  the  parameter  A is  defined  as 


A = 


^irS^2m 

h 


The  dependence  on  the  interface  potential  is  through  0,  the  average  of  the  two  inter- 
face potentials.  The  lack  of  dependence  on  A0,  the  difference  in  interface  potentials,  is 
reflected  by  the  lack  of  a term.  Apparently  this  lack  of  polarity  is  caused  by  the 
choice  of  expansion  about  an  average  Fermi  level.  Stratton  [7]  has  chosen  instead  to 
expand  about  the  Fermi  level  on  one  side.  Brinkman  et  al.  [8]  give  a low-voltage  expan- 
sion for  this  case  as 


where 


J = 


— 1/2 

me  .A<!> 
e 


2ir^h^ 


^eV~  ^ (eV)2  + 
^ 120 


_l/2 

160 


(eV)^ 


0 = 00  + 


Here,  not  only  is  the  cubic  term  3 times  larger  than  in  Eq.  (4),  but  there  is  a significant 
t quadratic  term  proportional  to  A0.  Unfortunately  this  expression  is  only  approximate, 
and  Brinkman  et  al.  do  not  state  the  approximations  involved.  Thus  the  simplest  type  of 
theory  is  somewhat  uncertain  even  for  the  second-  and  third-order  terms. 


JUNCTION  FABRICATION 

The  laboratory  fabrication  of  an  electron  tunneling  junction  is  Ulustrated  in  Fig.  3. 
The  procedure  involves  the  vacuum  evaporation  of  an  aluminum  strip  !^150  nm  thick  on- 
to a clean  glass  substrate.  The  strip  is  then  allowed  to  oxidize  at  room  temperature  from 
12  to  24  hours.  Following  this  oxidation  a second  similar  strip  is  deposited  at  ri^t 
angles  to  the  first  strip.  The  resulting  junction  sandwich  area  consist  of  two  parallel 
plates  of  thin  aluminum  separated  by  an  aluminum  oxide  insulating  film  approximately 
3 nm  thick.  Such  a structure  serves  as  a useful  tool  in  the  experimental  study  of  the 
electron  tunneling  phenomena  described  in  the  previous  section. 

In  practice,  five  to  ten  such  junctions  were  usually  fabricated  simultaneously  under 
the  same  conditions  to  establish  controls  and  statistics.  The  vacuum-chamber  arrange- 
ment for  evaporation,  masking,  and  monitoring  is  shown  in  Fig.  4.  Prior  to  evaporation 
a Vac-Ion  pump  is  used  to  evacuate  the  chamber  to  a pressure  of  5 X lO-"^  torr.  After 
degassing  the  tungsten  basket  filament  and  its  aluminum  contents  (0.5  g of  99.999+% 
aluminum),  the  aluminum  is  evaporated  rapidly  in  ^ 1.0  minute.  During  the  evaporation 
the  chamber  pressure  typically  increases  to  ^ 8 X 10"®  torr,  A shutter  arrangement  is 
used  to  control  the  starting  and  stopping  of  aluminum  deposition.  (The  circular  shutter 
is  displaced  to  the  left  in  Fig.  4.)  The  film  thickness  is  monitored  during  evaporation 
both  by  measuring  the  film  resistance  and  by  a direct-reading  digital  monitor  of  the 
crystal  oscillator  type  (Sloan  200).  The  programmable  monitor  is  also  used  to  follow 
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Fig.  3 — Laboratory  technique  used  in  fabricating  tunneling- 
junction  sandwich  structures  of  AI-AI2O3-AI.  The  aluminum 
oxide  film  is  grown  at  room  temperature  on  a vacuum-evaporated 
aluminum  strip  and  overlayed  with  a second  aluminum  strip  as 
shown. 


directly  the  subsequent  growth  of  the  aluminum  oxide.  After  the  aluminum  film  evapora- 
tion, the  monitor  is  reset  to  zero  thickness,  and  the  density  of  AI2O3  is  dialed  into  the 
monitor.  After  the  pressure  in  the  chamber  is  increased  to  1 atmosphere  of  pure  oxygen, 
the  AI2O3  thickness  increases  from  zero  to  ;^2  nm  in  the  first  5 minutes  and  then  grows 
more  slowly,  reaching  2.2  to  2.7  nm  after  about  1 hour.  Although  no  oxide  growth  was 
ever  observed  in  this  study  after  a few  hours,  the  second  evaporation  of  aluminum  over 
the  oxide  growth  was  not  carried  out  until  12  to  24  hours  later. 

Figure  5 shows  a typical  assembly  of  five  junctions  at  three  stages  of  fabrication. 
Initially  thin  strips  of  indium  are  “soldered”  onto  the  glass  microscope  slide  to  form 
electrical  contacts.  (Molten  indium  readily  adheres  to  a clean  glass  surface.)  A mask  is 
carefully  aligned  over  the  contacts,  and  a strip  of  aluminum  is  evaporated  onto  the  glass 
and  indium  contacts.  After  the  strip  is  oxidized,  a set  of  five  aluminum-film  strips  are 
similarly  evaporated  simultaneously  at  right  angles  to  the  first  strip.  Finally  No.  36 
copper  magnet-wire  leads  are  soldered  onto  the  indium  contacts. 


CURRENT-VOLTAGE  CHARACTERISTICS 

The  electrical  properties  of  these  junctions  which  are  of  primary  interest  in  this 
application  are  their  nonlinear  resistance  characteristics  and  their  generation  of  RF  inter- 
modulation.  The  nonlinear  DC  current-voltage  (I-V)  characteristics  of  such  junctions  were 
measured  with  a standard  four-point  arrangement  (Fig.  6).  A programmable  constant- 
current  supply  (Keithley  227)  is  used  to  drive  from  3 nA  to  1 mA  through  two  legs  of 
the  junction  film  leads.  A 5-1 /2-digit  digital  voltmeter  is  used  to  measure  the  voltage 
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Fig.  4 — Vacuum-chamber  arrangement  for  evaporation,  masking,  and  monitoring.  The  crystal  head  of  the 
Sloan  thickness  monitor  is  normally  mounted  in  the  plane  of  the  substrate  surfaces  and  is  not  shown  in  this 
photograph. 


developed  across  the  junction.  Since  the  input  impedance  of  the  voltmeter  is  >10®  ohms, 
negligible  current  from  the  current  source  flows  throu^  the  voltmeter  circuit,  and  the 
effects  of  lead  resistance  and  contact  resistance  are  essentially  eliminated;  that  is,  the 
only  IR  drop  seen  by  the  voltmeter  is  across  the  junction  area.  The  current  is  measured 
by  a digital  multimeter  (Keithley  164). 

Figure  7 shows  the  current-voltage  characteristic  curve  for  one  of  the  early  0.0168- 
cm^  devices.  On  the  scale  shown  the  experimental  curve  becomes  visibly  nonlinear  at 
about  80  mV  and  exhibits  a behavior  typical  of  all  the  junctions  measured.  On  the  basis 
of  tunneling  theory,  one  would  expect  the  curve  to  be  described  by  the  cubic  dependence 
of  the  form  shown  in  the  fourth-order  expansion  in  Eq.  (4).  A computer  program, 
FORLAN,  was  written  to  facilitate  rapid  calculation  of  the  tunneling  current  as  a function 
of  the  junction  parameters.  This  program  uses  the  complete  analytical  Forlani-Minnaja 
expression  together  with  the  near-linear  (low-voltage)  experimental  values  of  /,  V,  and  Gq 
to  determine  the  values  of  S for  selected  values  of  (^._The  program  can  then  generate 
full-range  characteristic  curves  for  different  values  of  so  as  to  obtain  a best  fit  to  the 
experimental  I-V  curve. 
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Fig.  5 — A typical  assembly  of  Five  junctions  shown  at  three  stages  of  construction:  (A)  Pattern  of  indium 
contacts  are  attached  to  clean  glass  substrate.  (B)  Cross  strip  of  aluminum  is  vacuum  deposited  and  oxidized 
in  pure  oxygen  at  room  temperature  for  !^24  hours.  (C)  Finally,  five  parallel  strips  of  aluminum  are 
vacuum  deposited  at  right  angles  to  the  cross  strip  to  form  five  junctions. 
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Fig.  6 — Standard  four-point  arrangement  for  measuring  the  current-voltage  i 

characteristics  of  the  electron  tunneling  junctions  • 
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DEVICE  NO.  7-31-S  NO.  1.  JUNCTION  1 AT  90  HRS 
AREA-=0.01flecm^ 

Ro*2.22K  OHMS 
C -0.0329 

r,- 2.3  OHMS  ,,C 


FROM  FORLANI-MINNAJA 
TUNNELING  EXPANSION  > 
V-0.86eV 
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l-GoV  + aGoV» 

^1  Go=4.5xlo^  mhot^V. 
f~  ir- 11.025  »olf»~* 


0 100  200  300 

V (mV) 

Fig.  7 — The  current-voltage  characteristic  curve  for  a 0.0168-cm2 
junction  measured  90  hours  after  fabrication.  The  equivalent  cir- 
cuit of  the  junction  is  shown  in  the  upper  right. 


The  quality  of  fit  to  the  experimental  data  was  found  to  vary,  depending  on  the 
particular  set  of  junctions.  In  some  cases  an  excellent  fit  could  be  obtained,  and  in 
other  cases  the  calculated  curvature  deviated  substantially  from  the  experimental  data 
beyond  100  mV. 

Using  the  best-fit  values  of  and  0,  one  can  also  calculate  the  tunneling  current 
from  the  expansion  in  Eq.  (4).  The  dashed  line  in  Fig.  7 shows  an  example  of  such  a 
calculation.  It  was  found,  however,  that  the  current-voltage  characteristics  of  all  of  the 
junctions  could  be  fitted  considerably  better  by  using  an  empirical  cubic  equation  of  the 
form  of  Eq.  (4):  I = GqV  + aGgV^.  Here,  when  experimental  values  of  Gq  and  / and  V 
at  100  mV  and  200  mV  were  used,  an  average  value  of  a could  be  obtained  to  give  a 
satisfactory  fit  to  the  data.  The  solid  line  in  Fig.  7 shows  an  example  of  this  procedure. 
The  a’s  calculated  from  the  tunneling  theory  in  Eq.  (4)  were  generally  20%  to  50%  lower 
than  the  empirically  determined  a’s.  Since  a is  a measure  of  the  nonlinearity  and  will 
later  occur  explicitly  in  the  IM  expression,  we  will  use  the  empirical  value  of  a as 
determined  for  each  junction. 
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VOLTAGE  (mV) 


Fig.  8 — The  behavior  of  junction  resistance  as  a function  of  voltage  and 
time.  The  solid  curves  are  generated  by  the  FORLAN  program  and 
reprint  a best  fit  to  the  experimental  points  for  the  indicated  values 
of  0 and  S. 


It  is  also  important  to  emphasize  that  the  junction  parameters  Gq,  a,  0,  and  S do 
not  remain  constimt  in  time.  As  indicated,  the  characteristics  shown  in  Fig.  7 were 
measured  at  90  hours  after  fabrication.  The  junction  resistance  as  a function  of  applied 
voltage  for  this  same  junction  is  shown  in  Fig.  8 with  time  as  a parameter.  Here  die 
solid  curves  are  a best  fit  to  the  experimental  points  usi^  the  FORLAN  program.  The 
curves  show  the  sensitivity  of  the  junction  resistance  to  0 and  particularly  to  the  insulator 
thickness  S. 

Although  the  marked  change  of  resistance  with  time  does  not  prevent  the  device 
from  being  well  characterized,  it  does  necessitate  the  additional  inconvenience  of  keeping 
a time  history  of  the  individual  junctions  and  of  measuring  the  intermodulation  within  a 
few  hours  of  any  characteristic  measurement.  A time  history  of  four  typical  junctions  is 
shown  in  Fig.  9.  These  four  junctions  were  fabricated  simultaneously  under  fte  same 
conditions.  The  low-voltage  resistance  Rq  of  each  junction  was  tracked  for  2000  hours 
or  more.  The  intermodulation  was  generally  measured  from  6 to  48  hours  after  comple- 
tion of  fabrication.  This  time  lapse  was  found  to  be  convenient  both  because  of  the  slow 
rate  of  resistance  change  and  because  the  initial  resistance  (300  f2  to  3 kA)  was  mote 
suitable  for  IM  measurements. 

An  additional  type  of  behavior  that  is  characteristic  of  the  junctions  can  be  seen  for 
two  of  the  junctions,  which  self-shorted  after  about  1300  hours.  Most  of  the  junctioits 
that  were  visibly  free  of  fabrication  defects  lasted  beyond  a few  hundred  hours.  In  a few 
cases  such  shorts  were  observed  to  repair  themselves  (return  to  a resistance  along  the 
expected  curve)  following  a voltage  spike  due  to  improper  switching  procedure.  Ihis 
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Fig.  9 — The  time  history  of  junction  resistance  Rq  for  four  typical  junctions. 
The  rate  of  resistance  change  increases  continuoudy  for  several  hundred  to 
several  thousand  hours,  followed  usually  by  an  abrupt  self-aborting  drop  to 
fractions  of  an  ohm. 


behavior  suggests  the  burnout  of  a localized  filamentary  short.  'Rie  mechanism  of  the 
time^eipendent  behavior  is  not  understood,  but  it  appears  to  be  due  to  a combination  of 
diffusion  at  the  interface,  contaminant  absorption,  ^ge  effects,  and  mechanical  rupture 
[9] . However,  all  of  the  observed  nonlinear  and  time-dependent  behavior  is  assumed  to 
realistically  simulate  the  numerous  microscopic  tunneling  points  of  real  metal-to-metal 
contacts. 


CAPACITANCE  MEASUREMENTS 

In  addition  to  the  nonlinear  resistive  characteristics  of  such  jtmctions  an  important 
property  which  dominates  the  RF  conduction  is  the  junction  capacitance.  Although  most 
of  the  junctions  fabricated  were  of  small  area  (;^  0.015  cm^),  their  ciqsacitance  is  relatively 
large  bMause  of  the  small  plate  separation  of  nm.  The  capacitances  of  the  junctions 
were  measured  at  1.0  MHz  with  a Boonton  direct  capacitance  bridge  (Model  75D). 

Since  the  thin-film  device  is  not  purely  capacitive,  the  bridge  will  see  the  equivalent 
circuit  of  the  device  (Fig.  7).  At  balance  the  bridge  and  device  admittances  Yr  and  Yp 
respectively,  or  reciprocally  the  impedances,  will  be  matched,  so  that  Yp  = Yp  and 

Yr^Gq  +jU}Cg.  (5) 

Here  the  conductance  Gp  and  capacitance  Cp  are  read  directly  from  the  bridge.  For  the 
typical  junction  parameters  shown  in  Fig.  7 and  for  frequencies  F > 1 MHz,  Rq  » Xp, 
v^ere  Xp  = l/2sfC;  therefore 
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‘D 


+ J 


Cj2c2 
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2c2 


(6) 


By  equating  the  real  and  imaginary  parts  of  Eqs.  (5)  and  (6),  one  can  obtain  the  device 
capacitance  C and  series  resistance  r,  in  terms  of  the  bridge  values  Gg  and  Cb. 


In  using  the  value  of  capacitance  thus  determined  one  is  assuming  that  the  value  of 
C is  not  a function  of  frequency  or  applied  voltage.  The  value  of  C showed  no  measurable 
change  with  applied  bias  voltage  up  to  300  mV.  Cursory  experiments  at  audio  frequencies 
also  showed  no  frequency  dependence  of  C.  The  use  of  a fixed  value  of  C for  the  junction 
implies  that  the  dielectric  constant  of  the  AI2O3  is  not  a function  of  voltage  and  that  the 
capacitive  part  of  the  junction  will  conduct  UHF  linearly.  A measurement  of  the  capacit- 
ance also  permits  a determination  of  the  effective  plate  separation  of  the  junction  sand- 
wich. For  a parallel-plate  capacitor  the  separation  S in  nanometers  is  simply 

KA 

S = 0.8854  ^ , (7) 

where  K is  the  dielectric  constant  for  AI2O3,  A is  the  plate  area  in  cm2,  c is  the 
capacitance  in  (iF.  For  the  junction  values  shown  in  Fig.  1,  A = 0.0168  cm2  q = 
0.0329  pF.  For  a value  of  K * 10.44  [10] , Eq.  (7)  gives  a plate  separation  of  4.7  nm. 

This  separation  is  considerably  larger  than  either  the  2.56  nm  resulting  from  a parameter 
fit  to  the  I-V  characteristics  (Forlani-Minnaja  expansion  in  Fig.  7)  or  the  2.25  nm  resulting 
from  the  digital  thickness  monitor  measurements  discussed  in  the  previous  section.  This 
discrepancy  has  been  observed  by  other  investigators  [11]  and  is  believed  to  result  from  a 
different  averaging  process  inherent  in  the  different  measuring  techniques. 

When  IM  effects  are  calculated,  the  bridge  values  of  capacitance  are  used  rather  than 
the  values  calculated  from  tunneling  or  weight-gain  measurements,  since  presumably  the 
reactive  RF  current  and  voltages  will  behave  in  accordance  witii  the  effective  value  of  C 
measured  at  1.0  MHz.  It  will  also  be  seen  in  the  following  section  that  the  RF  voltage 
developed  across  the  resistive  part  of  the  junction  is  an  extremely  sensitive  function  of 
the  capacitive  reactance. 


INTERMODULATION  MEASUREMENTS 

A block  diagram  of  the  basic  IM  test  facility,  together  with  the  equivalent  circuit  of 
a tunneling  junction,  is  shown  in  Fig.  10.  Each  transmitter  power  output  can  be  varied 
independently  from  0 to  100  W.  The  maximum  third-order  (290-MHz)  intermodulation 
generated  by  the  test  facility  alone  is  ;6-140  dBm  at  50  dBm  (100  W)  total  power  output 
from  the  diplexer.  The  tunneling  junctions  under  test  are  placed  between  the  diplexer 
output  and  150  m of  RG-214  coaxial  cable.  The  cable  approximates  an  infinite  trans- 
mission line  and  ideal  termination  [1].  In  the  case  of  the  tunneling  junction  devices,  the 
power  levels  were  adjusted  so  that  = F2,  and  the  total  power  input  Pj  to  the  junction 
was  usually  set  to  1.0  W.  The  total  power  input  of  any  device  was  kept  below  about  4 W 
because  of  the  limited  ability  of  the  thin-film  conduction  strips  to  dissipate  power. 
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Fig.  10  — NRL  intermodulation  test  facility  (operated  by  the  NRL  Satellite  Communica- 
tions Branch).  The  equivalent  circuit  of  a tunneling  junction  is  shown  in  the  normal  test 
configuration. 


As  shown  in  Fig.  9,  the  junction  resistance  changes  with  time,  making  it  necessary 
to  measure  the  1-V  characteristics  just  prior  to  IM  measurements.  Capacitance  values  are 
also  measured  prior  to  such  IM  test  runs,  even  thou^  the  capacitance  remains  essentially 
constant  with  time.  The  resistance  change  with  time,  however,  affords  the  opportunity 
to  measure  the  intermodulation  as  a function  of  junction  resistance. 

The  assembly  of  junctions,  usually  in  a set  of  five  as  shown  in  Fig.  5,  is  mounted  in 
a special  chassis  box  (Fig.  11).  To  minimize  the  generation  of  spurious  inteimodulation, 
the  aluminum  chassis  box  contains  no  ferromagnetic  materials  and  is  fitted  with  low-IM 
type-N  UHF  connectors  [1] . The  IM  power  level  of  the  chassis  box  was  tested  using  a 
No.  14  bare  copper  shorting  bar  solder^  between  the  terminals.  These  tests  consistently 
showed  IM  levels  of  ;^-140  dBm  at  an  input  power  of  1.0  W.  One  must  also  make  sure 
that  the  thin-film  structure  associated  with  the  junction  assembly  does  not  produce 
intermodulation.  A measurement  of  this  effect  was  made  by  using  the  film  cross-strip  on 
a set  of  five  junctions  as  a shorting  bar.  At  a 1.0-W  input  the  IM  level  was  measured  at 
>140  dBm  1 5 dB.  A thin-film  strip,  fabricated  and  mounted  in  the  same  way  as  the 
junction  assemblies,  also  showed  the  same  results.  Therefore  any  intermodulation  above 
the  residual  level  of  -140  dBm  ± 5 dB  would  be  due  to  the  junction  and  not  caused  by 
the  chassis  box,  connectors,  or  device  structure  external  to  the  junction. 

Because  of  the  impedance  mismatch  of  the  chassis  box  to  the  50-ohin  coaxial  cable, 
45%  of  the  power  was  observed  to  be  reflected  when  using  a copper  shorting  bar.  When 
an  assembly  of  junction  devices  is  mounted  in  the  chassis  box,  the  transmitted  power 
drops  from  55%  to  as  low  as  30%,  depending  on  the  size  of  the  junctions.  In  most  cases 
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Fig.  11  — Junction  assembly  mounted  in  aluminum  chassis  box.  Nonfer- 
romagnetic materials  and  low-IM  type-N  UHF  connectors  are  employed  to 
minimize  spurious  generation  of  intermodulation. 


the  impedance  of  the  junction  and  leads  is  less  than  20  n.  The  transmitted  power  Pj  as 
experimentally  measured  is  used  in  calculating  the  RF  current  passing  throu^  the  junc- 
tions. 

To  calculate  the  IM  power  to  be  expected  from  a tunneling  junction,  we  return 
to  Eq.  (2).  For  the  tunneling  junctions  studied,  it  was  found  that  the  experimental  DC 
current-voltage  characteristics  can  be  adequately  described  by 

I = GqV^G^V^.  (8) 

In  thi-s  case  there  is  no  term,  since  the  tunneling  current  is  essentially  symmetric 
about  the  origin  of  the  /-V  characteristics.  Also,  as  can  be  seen  from  Eq.  (3),  any  V 
term  would  produce  no  IM  frequencies  falling  within  the  receive  band.  Equation  (8)  can 
be  written  in  a more  convenient  form  by  letting  a = G2/G0,  so  that 

I^GqV^oiGqV^,  (9) 

where  Gq  is  the  usual  linear  conductance  and  a is  a ratio  that  measures  the  departure 
from  linearity.  In  the  case  of  an  RF  input  voltage  of  the  form  of  Eq.  (1),  the  current  in 
Eq.  (9)  can  be  written  as  shown  in  Eq.  (3),  except  that  here  only  the  experimentally 
observed  angular  frequencies  tOj,  0)2  and  2oj2  - are  retained: 
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^(0  - Gq(Fj  cos  Wjf  + V2  COS  W2O  + - aGfjVi  V|  cos  (2cj2  - coj  )< 


(10) 


It  is  impoi^t  to  note  that  I(t)  is  the  RF  current  throu^  the  nonlinear  resistance  Rn 
thaWp(t)  IS  the  primary  current  at  frequencies  and  A,  and  that  /,w(f)  U the  inter- 
modulation  current  at  frequency  2/^2  - /j.  The  peak  IM  current  is  just 

g 

^IM(peak)  ” 4 (11) 


TTie  value  of  Vj  and  1^2  can  be  calculated  from  the  total  transmitted  povirer  /V-  The 
t^Emitted  power  Pj,  results  from  both  power  inputs  Pj  and  Pg  running  simultaneously 
at  frequencies  /j  and  ^2-  The  circuit  parameters  of  the  equivalent  circuits  in  Figs.  7 and 
10  have  the  following  ranges,  depending  on  junction  area  and  history: 

300n<P(j  <100k^2, 

0.026  n<Xc<  0.296  12  at  290  MHz, 

2 J2  < < 25  J2, 

Tj  = 50  J2. 

For  all  cases  Rq  » Xg,  so  that  essentially  all  of  the  transmitted  current  /j.  (fj'(t)  = 

'1  + /2  cos  0)20  passes  throu^  the  capacitance  C and  the  total  transmitted 

l^wer  Pj.  IS  dissipated  in  the  series  resistance  r*  and  the  load  resistance  r^.  Since  r,  + r. 
determines  the  current,  the  measured  RF  power  Pj.  is 

Pp  = < /f.(f ) > (r,  + Tg ),  (12) 

where 

•'0 

For  Pj  = P2  and  /j  = /j. 


(/j  cos  Wjt  +/2  cos  0)20^  dt. 
Eq.  (12)  yields 

/ 

T(peak)  L + 


Since  essentially  all  of  this  current  passes  throu^  the  capacitance,  the  voltages  V,  and  Vo 
developed  across  the  capacitance  C and  also  across  the  nonlinear  resistance  Rn  are 
approximately  '' 
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Vi=^Cl 


where  ^C1  and  Xq2  8re  the  capacitive  reactances  at  and  ^2  • these  values  of 

and  V2,  Eq.  (11)  becomes 


AM(peak) 


■ T/2 

rrr 


(13) 


This  IM  current  is  viewed  as  being  generated  in  the  nonlinear  element  Rq,  and,  as  can  be 
seen  from  Fig.  10,  this  current  source  is  paralleled  by  the  junction  capacitance  C and  the 
external  load  r,  + Tj.  Again  X^,  « r,  + at  290  MHz,  so  that  most  of  this  current  is 
shunted  throu^  C,  and  the  voltage  developed  across  C and  also  across  the  external  load 
r,  + fg  is 


'^IM(pcak)  ^IM(peak)-^C.IM’ 


(14) 


where  X^  m is  the  capacitive  reactance  at  290  MHz.  The  IM  power  dissipated  in 
is  then 


^m(pcak) 


(2c02  - Wi)t  > 


1 V 

1.  *^IM(peak) 

2 r.  + r„ 


(15) 


Combining  Eqs.  (13),  (14),  and  (15),  we  obtain  finally 


^IM  ~ 32 


y*  y4  y2 

^C1^C2^C,1M 


(16) 


This  equation  expresses  the  IM  power  level  in  terms  of  readily  measurable  junction  and 
test  parameters.  The  RF  test  facility  is  set  up  to  measure  the  IM  power  level  in  dBm 
units  (decibels  referenced  to  1 mW),  so  that  by  definition  Eq.  (16)  can  be  expressed  as 


PiM  (watts) 

(dBm)j^  = 10  log  -4 • (17) 

10"®  (watts) 


This  dBm  level  of  intermodulation  is  measured  directly  at  the  spectrum  analyzer  by 
comparison  to  a reference  signal,  as  shown  in  Fig.  10. 

The  junctions  initially  measured  in  the  IM  test  facility  had  relatively  large  areas 
{A  ^ 0.015  cm^)  and  were  operated  at  an  input  power  of  % 1.0  W.  For  a transmitter 
power  Pt<  * 0.5  W these  junctions  showed  intermodulation  levels  of  -110  to  -135  dBm, 
depending  on  the  values  of  Rq  and  C of  the  particular  junction.  Figure  7 shows  the  I-V 
characteristics  and  parameters  typical  of  this  type  of  junction.  If  one  uses  the  values 
shown  in  Fig.  7 together  with  Pj  » 0.5  W and  R^  = 50  ^2,  a calculation  of  the  IM  power 
from  Eq.  (16)  gives  a much  lower  IM  level  of  -239  dBm. 
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In  the  preceding  calculation  it  was  assumed  that  the  RF  current  is  uniformly  con- 
ducting through  the  film  strip  and  junction  area.  At  input  frequencies  of  250  and  270 
MHz,  however,  one  mi^t  expect  UHF  skin  effects  to  significantly  alter  the  IM  level.  An 
accurate  calculation  of  the  RF  current  distribution  and  resultant  intermodulation  is  an 
extremely  complex  problem  for  the  geometry  of  the  thin-film  junction  device. 

One  simplifying  approximation,  which  has  been  employed  by  Chapman  et  al.[2] 
regarding  waveguide  flanges,  is  to  restrict  the  RF  current  across  the  tunneling  area  to  the 
same  depth  as  the  skin  depth  in  the  adjacent  conductor.  Although  the  geometry  of  a 
thin-film  device  differs  significantly  from  a comparatively  large  flange,  a similar  restriction 
arises  through  edge  concentration  of  RF  current  in  thin-ribbon  conductors.  (Evidence  of 
such  a current  distribution  was  indicated  in  one  set  of  junctions  in  which  the  RF  power 
had  exceeded  10  W.  Small  thermally  produced  blisters  showed  a marked  concentration 
along  the  outer  edge  of  the  conducting  strips.) 

For  this  approximation  the  effective  RF  current  is  restricted  to  an  average  skin 
depth  of  5.16  X 10"“*  cm  around  the  perimeter  of  the  actual  junction  area.  For  the 
junction  parameters  shown  in  Fig.  7,  the  effective  RP  conduction  area  is  2.67  X 10"'^ 
cm^,  or  about  1/63  the  apparent  area.  Thus  the  effective  junction  resistance  and  capaci- 
tance become  ^^1.39  X 10^  SI  and  ;^5.23  X 10“^^  P respectively.  With  these  effective 
values  of  Rq  and  C,  Eq.  (16)  gives  an  intermodulation  level  of  -131  dBm,  which  is  within 
the  experimental  range  observed. 

To  raise  the  IM  level  and  to  examine  any  area  effects,  several  sets  of  five-junction 
assemblies  were  fabricated  in  graded  sizes  ranging  from  0.015  cm^  down  to  0.0015  cm^. 

It  was  found  that  the  size  scaling  effect  on  capacitance  C and  junction  resistance  Rq 
essentially  followed  the  expected  variation  for  such  junction  elements  in  parallel.  As  the 
junction  area  is  decreased,  one  would  expect,  according  to  Eq.  (16),  the  capacitive  react- 
ance product  (Xci^c2^c  IM  ^ increase  much  faster  than  Rq  and  to  result  in  a 

rapid  increase  in  with’  decreasing  area.  This  increase  is  largely  offset,  however,  by  an 
increase  in  the  ratio  of  effective  RP  conduction  area  to  total  junction  area  as  the  junction 
area  decreases.  When  the  present  fabrication  techniques  are  used,  a reduction  in  junction 
area  by  more  than  an  order  of  magnitude  becomes  impractical,  and  the  power-handling 
abUity  becomes  too  low  to  conveniently  make  IM  measurements. 

Figure  12  shows  typical  results  of  IM  measurements  on  a set  of  five  junctions  ranging 
in  size  from  0.0143  cm^  to  0.0018  cm2.  The  measured  value  of  the  junction  capacitance, 
is  shown  along  the  abscissa.  The  IM  power  level  in  dBm  units  and  in  watts  is  shown  on 
the  left  and  right  ordinates  respectively.  The  input  power  Pj  is  held  constant  at  +30  dBm 
(1.0  W).  Over  the  capacitance  range  examined  the  transmitted  power  Py  varies  from 
0.39  W to  0.51  W but  remains  relatively  constant  on  ttie  scale  shown.  The  IM  level  for 
each  junction  is  plotted  as  a square  point  along  with  the  estimated  error  bars  of  ±3  dB. 
These  experimental  points  are  connected  with  a smooth  dashed  curve.  The  solid  curve 
shows  the  intermodulation  calculated  from  Eq.  (16)  using  the  experimentally  measured 
junction  parameters  with  the  values  of  Eg  ^ corrected  to  the  effective  RF  skin  depths. 
The  uncorrected  junction  parameters  are  shown  in  Fig.  13.  Similar  sets  of  junctions 
showed  the  same  functional  dependence  of  intermodulation  on  area  and  capacitance  as 
illustrated  in  Fig.  12. 
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Fig.  12  — Experimental  and  calculated  intern- adulation 
power  levels  generated  by  electron  tunneling  as  a function 
of  junction  capacitance.  Here  the  variation  in  capacitance 
results  from  a gradation  in  junction  areas  ranging  from 
0.0018  cm^  to  0.0143  cm^.  All  junctions  were  fabricated 
simultaneously  and  have  the  same  oxide  thickness  of  ^2.5 
nm.  The  junction  parameters  and  /-V  characteristics  for 
each  of  the  Rve  junctions  are  shown  in  Fig.  13. 


Althou^  the  intermodulation  calculated  in  the  preceding  manner  followed  the 
exjperimental  variation  with  junction  parameters  fairly  well,  the  absolute  magnitude  of 
the  calculated  IM  curve  was  found  to  be  shifted  from  the  measured  intermodulation  by 
±15  dB.  This  agreement  is  perhaps  as  good  as  can  be  expected  in  view  of  the  crude  RF 
skin-depth  approximations  used.  Additional  smaller  effects  which  are  not  taken  into 
account  by  Eq.  (16)  are  distributed  resistance  and  capacitance  effects,  nonlinear  capaci- 
tive effects,  and  frequency-dependent  effects. 

The  equivalent-circuit  analysis  leading  to  Eq.  (16)  tacitly  assumes  that  the  film 
resistance  over  the  junction  area  is  zero.  Actually  this  resistance  is  generally  larger  than 
the  capacitive  reactance  at  270  MHz;  for  example,  in  junction  1 of  Fig.  13,  Xq  = 

0.0286  £2  and  tf  0.185  £2.  Thus,  one  should  consider  the  effects  of  distributed  resis- 
tance and  capacitance  along  the  junction  area.  With  use  of  a two-dimensional  model,  both 
an  exact  continuous-distribution  analysis  and  a numerical  discrete  analysis  showed  a 
voltage  distribution  V(x)  that  varied  significantly  along  the  width  of  the  junction.  In 
addition  the  RF  edge  conduction  along  the  film  strip  leading  up  to  the  junction  area  will 
increase  the  effective  value  of  r, . A calculation  of  ftese  two  combined  effects  for  junc- 
tion 1 of  Fig.  13  using  a discrete  four-segment  equivalent  circuit  yielded  a value  of  ^IM 
which  differed  from  the  previously  corrected  value  by  only  1 dB.  However,  this  small  net 
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Fig.  13  — Current-voltage  characteristics  as  a function  of  area 
and  capacitance  for  the  five  junctions  measured  in  Fig.  12. 
The  corresponding  junction  parameters  are  shown  in  the  upper- 
left  table. 


correction  is  due  to  the  fortuitous  approximate  cancellation  of  the  two  effects  for  the 
geometry  chosen. 

llie  a of  Eq.  (16)  is  derived  from  the  nonlinear  DC  resistance  characteristics  of  the 
junction.  Hiis  nonlinearity  is  presumed  to  arise  solely  from  electron  tunneling,  as 
described  previously.  However,  known  nonlinear  capacitance  effects  can  arise  through 
nonlinear  dielectric  and  electrostriction  effects  [12].  Electrostriction  can  also  change  the 
tunneling  resistance  throu^  a volume  change.  A rou^  calculation  of  the  volume-change 
effect  predicted  an  intermodulation  well  below  the  level  of  detectability. 

Equation  (16)  is  also  applied  to  RF  power  measurements  of  and  Pf  in  the  UHF 
range  (240  MHz  to  290  MHz),  whereas  the  parameters  a,  Rq,  and  r,  are  determined  from 
DC  measurements  and  C is  measured  at  1 MHz.  The  assumption  of  frequency 
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independence  seems  reasonable,  since  the  resistivity  and  dielectric  constant  of  AlgOg 
appear  to  be  relatively  constant  over  a wide  range  of  frequencies  [13,14]. 


ION  IMPLANTATION  EFFECTS 

To  eliminate  or  reduce  the  IM  effects  of  tunneling  junctions  in  metal-to-metal  con- 
tacts, one  apparently  must  either  eliminate  the  surface  insulating  films  (oxides  or  other- 
wise) or  basically  modify  the  structure  of  the  metal  oxide  surface.  Fabrication  of  full- 
scale  metal  hardware  that  is  essentially  free  of  metal  oxide  films  both  at  the  surfaces  and 
under  metal  platings  such  as  gold  would  appear  to  offer  a reasonable  solution.  Alterna- 
tively, modification  of  the  metal  oxide  surfaces  should  ideally  make  any  contacts  inore 
conductive,  linear,  and  free  of  time-dependent  effects.  We  have  explored  implantation  of 
metallic  ions  in  the  oxide  surface  as  a possible  way  to  achieve  the  desired  results. 

Several  tunneling  jxmctions  were  fabricated  similarly  as  previously  described,  except 

ions  of  2 keV  energy  were  implanted  into  the  AI2O3  surface  prior  to  the  evapora- 
tion of  the  second  overlapping  aluminum  strip.  Theoretical  estimates  indicate  that  2-keV 
Ag'*’  ions  have  a mean  range  in  AI2O3  of  nm.  Ideally  these  silver  atoms  will  have  a 
Gaussian  distribution  about  the  mean  range  with  a standard  deviation  of  ^ 0.6  nm.  A 
plot  of  Gaussian  distribution  shows  range  straggling  extending  from  the  surface  to  ^ 3.2  nm, 
with  about  12%  of  the  silver  atoms  stopping  in  the  aluminum  beyond  the  AlgOg. 

The  actual  distribution  of  the  silver  atoms  in  the  AI2O3  and  A1  may  be  quite  dif- 
ferent than  just  described  because  of  three  effects  that  have  not  been  measured. 

• Uncertainty  in  the  low-energy  range  of  Ag"^  ions  in  AI2O3  (existing  data  indicate 
that  heavy-ion  ranges  in  light  substrates  may  be  as  much  as  twice  the  theoretical  estimates), 

• Removal  of  the  AI2O3  surface  by  ion  sputtering,  and 

• Oxide  regrowth  when  the  film  is  removed  from  the  vacuum  environment. 

To  evaluate  the  electrical  characteristics  of  such  implanted  junctions,  standard 
reference  monitor  junctions  were  also  fabricated  simultaneously  on  the  same  substrate 
and  under  the  same  conditions  without  ion  implantation.  All  electrical  measurements  are 
then  referenced  to  the  monitor  junctions.  This  procedure  is  necessary  because  different 
batches  of  junctions  will  generally  have  different  values  of  Rg  because  the  electrical 
characteristics  will  be  measured  at  different  times. 

Initially  two  sets  of  junctions  (three  reference  junctions  per  set)  were  implanted. 

One  set  was  implanted  with  a fluence  of  atoms/cm2,  and  the  other  set  was 

implanted  with  a hi^er  fluence  of  atoms/cm2.  The  junction  implanted  at  the 

hi^er  fluence  showed  a sharp  drop  in  resistance  from  an  average  resistance  of  36.3  kS2  to 
an  average  resistance  of  13.5  or  a ratio  of  implanted  resistance  RoiAg)  to  reference 
resistance  R-(Ref.)  of  ^4  X 10-4.  -nie  junctions  implanted  at  lower  fluence,  however, 
unexpecteiy  showed  an  increase  from  12.3  kfl  to  23  Mf2,  or  a ratio  JiQ(Ag)/i!q(Ref.)  * 
1.9  X 10®.  The  different  implant  fluences  produced  opposite  effects  of  change  in 
resistance  that  differed  by  over  six  orders  of  magnitude! 
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Fig.  14  — Effect  on  junction  resistance  of  implantation  of  2-keV 
silver  ions  in  the  AI2O3  surface  as  a function  of  fluence  of  silver 
atoms.  The  ratio  of  ion-implanted  resistance  Rg(Ag)  to  an  identical 
unimplanted  resistance  Rg(Ref.)  is  plotted  on  the  ordinate. 


In  view  of  this  surprising  result  a more  elaborate  set  of  junctions  were  fabricated  to 
facilitate  measurement  of  the  junction  resistance  as  a function  of  the  fluence  of  silver 
atoms  implanted.  Figure  14  ^ows  the  results  of  these  measurements.  The  ratio 
i2Q(Ag)/i2o(Ref.)  is  plotted  on  the  ordinate,  and  the  fluence  of  silver  atoms  and  of 
charge  are  plotted  along  the  top  and  bottom  abscissa  respectively.  The  solid  line  con- 
nects points  measured  immediately  after  fabrication,  and  the  dashed  line  connects  points 
measured  about  2 hours  later.  The  resistances  are  changing  rapidly  during  this  initial 
period,  but  the  ratio  remains  relatively  constant  in  time.  The  two  square  points  are  from 
the  original  measurements  and  are  seen  to  be  consistent  with  the  more  detailed  later 
results.  The  error  bar  on  one  square  point  is  significantly  hi^er  than  the  errors  for  the 
other  points  and  resulted  from  an  equipment-produced  uncertainty  in  the  charge  fluence 
at  this  one  point.  The  increase  in  junction  resistance  for  implants  up  to  4 X 
atoms/cm^  was  completely  unexpected,  and  the  detailed  physical  mechanisms  responsible 
for  the  functional  behavior  of  /io(Ag)/i{Q(Ref.)  with  the  silver  implant  fluence  are  not 
yet  fully  understood.  Recent  work  [15]  shows  that  ion  implantation  in  the  oxide  sur- 
face of  metals  has  the  effect  of  increasing  the  oxide  thickness.  This  effect  would  account 
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for  the  initial  increase  in  tunneling  resistance.  As  the  silver  density  is  increased,  the  oxide 
conductivity  apparently  increases  faster  than  the  competing  thickness-growth,  resulting  in 
a low-resistance  junction.  For  those  junctions  implanted  to  fluences  of  > 10^®  atoms/cm^ 
the  resistance  fig  was  generally  reduced  to  a few  ohms  or  fractions  of  an  ohm.  The 
current-voltage  characteristics  of  the  implanted  junctions  showed  better  linear  behavior 
than  unimplanted  junctions  over  the  same  current  range.  Although  these  junctions  had 
greatly  improved  conduction  characteristics,  they  continued  to  show  the  familiar  time- 
dependent  increase  of  fig • 

Preliminary  test  of  these  junctions  at  the  IM  test  facility  confirmed  that  the  low- 
resistance  implanted  junctions  did  indeed  have  a low  intermodulation  of  -135  dBm. 

For  such  an  implantation  technique  to  be  useful  on  full  scale  hardware,  however,  one 
would  have  to  solve  the  formidable  problems  of  the  time-dependent  behavior  and 
mechanically  fragile  surfaces. 


CONCLUSIONS 

From  our  experimental  and  theoretical  investigation  of  the  generation  of  intermodula- 
tion by  electron  tunneling  in  AI-AI2O3-AI  junctions,  we  list  the  following  conclusions: 

• Tunneling  junctions  ranging  in  areas  from  0.015  cm^  down  to  ^0.0015  cm2 

showed  IM  levels  of  150  dB  to  ^ 110  dB  below  the  total  transmitted  power  of  0.5  W. 
This  corresponds  to  a PiulPm  respectively.  Since  ideally  the 

ratio  Fjuj/Pgig  should  be  10-19,  the  electron  tunneling  mechanism  is  seen  to  be  a 
significant  source  of  intermodulation. 

• The  functional  dependence  of  the  IM  power  on  the  effective  RF  junction  parame- 
ters appears  to  be  satisfactorily  described  by  Eq.  (16). 

• Although  the  main  featiures  of  the  current-voltJ^e  characteristics  of  the  junctions 
are  described  by  the  Forlani-Minnaja  theory,  an  accurate  fit  to  the  experimental  data  was 
not  always  possible  over  the  entire  range  of  experimental  values. 

• No  significant  effects  of  hysteresis  in  the  I-V  characteristics  were  observed.  How- 
ever, a small  resistance  asymmetry  was  always  observed  at  the  higher  voltages.  If  the 
film  strip  oxidized  first  was  operated  at  a positive  polarity,  the  junction  resistance  was 
lower  than  the  reverse  polarity  by  3%  to  7%.  TTius  the  junctions  exhibit  a small  rectify- 
ing component  that  is  thought  to  arise  from  the  asymmetry  of  the  AI2O3  density  profile 
at  the  two  interfaces. 

• In  well-characterized  AI-AI2O3-AI  junction  devices  the  time-dependent  behavior  is 
prominent,  with  the  junction  resistance  increasing  continuously  for  several  hundred  to 
several  thousand  hours,  followed  usually  by  an  abrupt  self-shorting  drop  to  firactions  of 
an  ohm.  During  such  resistance  increase,  the  capacitance  remains  essentially  constant. 
Such  time-dependent  behavior  is  assumed  to  take  place  randomly  in  the  numerous  micro- 
scopic tunneling  points  in  real  hardware  contacts.  This  is  a mechanism  that  could  contri- 
bute to  the  erratic  fluctuation  of  intermodulation  observed  in  macroscopic  contacts. 
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• Ion  implantation  of  ^ 10^®  silver  atoms  per  cm^  in  the  oxide  surfaces  produced 
low-resistance  and  low-IM  junctions  but  failed  to  stop  the  characteristic  increase  of  junc- 
tion resistance  with  time.  Our  initial  implantation  results  are  sufficiently  promising, 
however,  to  warrant  further  investigation  of  the  possible  beneficial  effects  of  employing 
different  metallic-ion  species,  greater  implantation  depths,  and  hi^er  fluences. 
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